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We investigate heat transport in a spin-i Heisenberg chain, coupled locally to independent thermal 
baths of different temperature. The analysis is carried out within the framework of the theory of open 
systems by means of appropriate quantum master equations. The standard microscopic derivation 
of the weak-coupling Lindblad equation in the secular approximation is considered, and shown 
to be inadequate for the description of stationary nonequilibrium properties like a non-vanishing 
energy current. Furthermore, we derive an alternative master equation that is capable to describe 
a stationary energy current and, at the same time, leads to a completely positive dynamical map. 
This paves the way for efficient numerical investigations of heat transport in larger systems based 
on Monte Carlo wave function techniques. 

PACS numbers: 05.60.Gg, 05.30.-d, 05.70.Ln 



I. INTRODUCTION 

The transport behavior of one-dimensional systems 
has intensively been investigated both in classical and 
in quantum mechanical context for several decades now. 
In the classical domain it seems to be largely accepted 
that normal energy transport, i.e., spatial diffusion in- 
stead of ballistic transport or localization requires the 
chaotic dynamics of a non-integrable system [l],[2| . In the 
non-classical regime the question whether normal trans- 
port behavior may arise from the underlying quantum 
mechanical equations of motion is a controversial issue 
@i Hi ■ This is mostly due to the nontrivial character of 
the question, how energy or heat is transported through 
a system on a microscopic level of description. 

Among the most prominent techniques of theoretical 
description is the use of the Green-Kubo Formula 0, 0] 
that is derived on the basis of linear response theory. 
Here, similar to the classical case, strong evidence arises 
that integrability leads to diverging transport coefficients 
and that thus normal transport is not to be expected 
[|. An overview of theoretical approaches along with 
their respective subtleties, pitfalls and shortcomings was 
recently given in Q. 

The article at hand addresses the problem within 
the framework of the theory of open quantum systems 
d, Qil [H| by coupling the system of interest explicitly to 
environments of different temperature. In recent years 
many investigations on heat tr ansp ort have been car- 
ried out in this framework [l^, O EH- According to 
the results of these studies numerical evidence for normal 
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transport behavior arises, even in the integrable system. 
However, the investigations are generally limited to small 
system sizes as the computation is highly demanding. 

To investigate a real bulk property of the material, 
e.g., the heat conductivity, one primary aim is to apply 
methods that can cope with larger systems in order to 
rule out that the mentioned findings are merely due to 
finite size effects. The present article constitutes one step 
towards this aim. 

The common setup in models of open quantum sys- 
tems comprises the system of interest and an environ- 
ment. Employing various approximation schemes for the 
system-environment coupling one derives effective dissi- 
pative equations of motion for the reduced density ma- 
trix ps of the open system. These equations are called 
Markovian quantum master equations (QME). Here, the 
term Markovian refers to the special time-local structure 
that arises when all memory effects are negligible. The 
QME is usually required to be in Lindblad form in order 
to guarantee the preservation of the normalization and 
of the positivity of the reduced density matrix, as well as 
the complete positivity of the resulting dynamical map 

EE El. 

Once equipped with a QME of Lindblad structure the 
way is paved for efficient numerical studies of larger sys- 
tems by means of the stochastic simulation techniques 
provided by the Monte Carlo wave function method (see, 
e.g., [1, [13] and references therein). 

In contrast to the common open system setup de- 
scribed above, where the environment serves as a heat 
bath driving the system into thermal equilibrium, the 
introduction of a second heat bath of different tempera- 
ture leads to dynamics that may feature a nonequilibrium 
steady state. Such states typically feature both energy 
currents as well as temperature gradients. Unfortunately, 
in the present case the derivation of a QME in Lind- 



2 



blad form, that reasonably describes the expected physi- 
cal behavior turns out to be other than straightforward. 
Lacking any QME in Lindblad form also the mentioned 
efficient stochastic simulation technique is no longer ap- 
plicable. As an example let us study energy transport in 
one-dimensional Heisenberg spin-i chains that are cou- 
pled locally to heat baths of different temperature. 

In order to illuminate the problem that arises in the 
attempt to apply the standard recipe to transport dy- 
namics we shall focus on one of the aforementioned ap- 
proximations that are invoked in the derivation of a QME 
which actually does feature Lindblad form. The most 
common approximation which yields a Lindblad form (in 
the weak- coupling limit) is the so-called secular approx- 
imation (SA) [18|, [l9|. Essentially, this approximation 
consists in replacing the generator of the interaction pic- 
ture master equation by its time average. The SA is 
justified as long as the time scale set by the differences 
of the Bohr frequencies of the open system is short com- 
pared to the relaxation time. In many cases it is equiva- 
lent to the rotating wave approximation, but there exist 
counterexamples |20j. 

We observe that, once the SA has been performed, one 
is left with a QME that predicts a vanishing energy cur- 
rent in the presence of a finite temperature gradient. The 
standard weak-coupling Lindblad master equation that is 
obtained on the basis of the SA is therefore inappropri- 
ate to describe heat transport. Hence, in order to be 
able to apply the standard Monte Carlo wave function 
techniques to heat transport one needs an appropriate 
quantum master equation which leads to a finite station- 
ary energy current and, at the same time, is of Lindblad 
form. Here, we derive such a master equation which is 
valid in the regime of weak internal couplings. 

The paper is organized as follows. In Sec.Qllwe intro- 
duce the details of our model and give the definition of 
the energy current operator. Moreover, we briefly recall 
the derivation of the Redfield master equation for the 
density matrix of the spin chain. Section IIIII contains a 
discussion of the SA and of its influence in the descrip- 
tion of heat transport. We further derive a completely 
positive Markovian master equation for the dynamics of 
the spin chain and demonstrate that it yields an excel- 
lent approximation of the dynamics given by the Redfield 
equation. Several conclusions and implications for future 
investigations are presented in Sec. HVl 



II. THEORETICAL BACKGROUND 

A. The Model and the current operator 

We investigate heat transport in spin-^ chains of 
length n, coupled to heat baths of different temperature 
at both ends (see Fig. [I]). The Hamiltonian of the spin 
chain reads 



Ha — H] 
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FIG. 1: Spin-i chain coupled to heat baths of different tem- 
perature. The stationary state features a non-vanishing en- 
ergy flux through the chain of spins from the hotter towards 
the colder heat bath. 



The first term describes a local energy contribution due 
to an external field of strength £1, 



(2) 



whereas the second term accounts for a homogeneous 
nearest-neighbor Heisenberg interaction with a coupling 
parameter A, 



V 



^2 t/(a*,m+i) 



(m+i) 



(3) 



& = ((J x ,a y ,a z ) T is the spin vector operator with the 
well-known Pauli spin matrices as its elements. Braced 
upper indices label the respective spin. 

In order to obtain an operator for the energy current 
between two adjacent spins in our system, we consider 
the time evolution of the local energy operator given by 
the Heisenberg equation of motion for operators at site 
M 



d H\oc _ -rfr rVGOi , 9 MiA 
-i[H s ,H loc ] + -H loc 



(4) 



Since is not explicitly time dependent Eq. (j4]) be- 

comes after inserting Eq. ([!]) 



00 



dt 



(\y^ 1 *\Hj£] + \y^i\Hg>]) . ( 5 ) 

Assuming that the local energy is a conserved quantity, 
which is justified when f2 is large compared to A, we can 
rewrite Eq. §5§ as 



dff, 



loc 



df 



divJ = - jO*- 1 ./*) 



(6) 



where we introduced a discretized version of the continu- 
ity equation. By comparing Eqs. ([5]) and ([6]) we deduce 

being the energy current flowing from site fi to site ji + 1. 
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B. Redfield Equation 

The total system, i.e., system and both heat baths is 
described by the Hamiltonian 



2 

-fftot = H s + 



Mi) , 6-W 



(8) 



iJg stands for the Hamiltonian of the ith bath, 
denotes the respective system-bath interaction. 

In the following we will first turn to a scenario with a 
single bath for the purpose of a compact notation. The 
time evolution of the total system is described by the von 
Neumann equation which reads 



d/3to 
dt 



— — i[Ht t, ptot] 



(9) 



Ptot represents the density operator of the total system. 
The aim is to derive an equation of motion for the density 
operator pg of the spin chain which is defined through the 
partial trace taken over the bath variables, 



Ps(t) = Tr B {Ptot(t)} 



(10) 



To this end, we invoke the Born-Markov approximation 
and assume that the bath is in a thermal equilibrium 
state at inverse temperature /3. 

If the interaction Hamiltonian is cast into the form 



(11) 



with Hcrmitian operators X k and Y k of system and bath 
respectively, the standard procedure (cf. @) yields the 
following Schrodinger picture QME 



^h = -i[H s ,p s ]+V(p s ). 



(12) 



The first term describes the free dynamics, whereas the 
dissipative part 



E 

kl 



&TT lk ( T , (3)[X k (-T) / S s ,i'j]+H.c. 



(13) 



accounts for the influence of the environment. The time 
dependence of the operator X k (— r) in Eq. ([13]) has to be 
interpreted as 



X k {-r) = e 



-iH S T 



X k e 



(14) 



The bath inverse temperature (3 is embodied in the bath 
correlation functions 



T h (t,0) = (Y k (r)Yi) B = Tr B {p B Y k (r)Yi} (15) 



with 



-PH B 



Pb 



(16) 



Equation ([12]) is known as Redfield master equation, and 
has many applications ranging from NMR to the descrip- 
tion of chemical dynamical systems [111, [^ ■ We will use 
this equation in the following as a reference point to judge 
whether further approximations are suitable or not. 



C. Local Environment Coupling 

The chain of n spins shall be coupled only locally to the 
respective baths, i.e., via the outermost spins (cf. Fig.[T]). 
Consider, e.g., the system operator that couples to the 
left heat bath 



x L = &^ ® i< 2 > ® • • • ® iw = 



(17) 



Thus, the double sum in Eq. (|13p contains only one single 
term. The corresponding Redfield dissipator reads 

/>oo 

PlGSs)=/ drr(r,/3 L )[X L (-T)p S) l L ]+H.c. (18) 
Jo 

We choose a simple bath modelled by an infinite number 
of independent harmonic oscillators, 



H B = ^2^j k b\b k , 



(19) 



k=l 



with the usual bosonic creation an annihilation operators 



bt, b k . The system-environment coupling as defined in 



Eq. (TiT]) is taken to be linear in the oscillator amplitudes, 

oo 

f L =J> fc 6t +c *i fcj (20) 



k=l 



where the c k are coupling constants. Plugging Eq. (|20|) 
into Eq. (fT5|) the bath correlation function is found to be 



r(r,/3 L ) 



&lu e iwT r(w,/3 L ) 



(21) 



with 



r(w, /3 L ) = - [J(w) - J(-w)] N(u, fa) . (22) 



Here, J{ui) denotes the spectral density, N(w,fa) repre- 
sents the Planck distribution 



1 



(23) 



and k is a coupling parameter. Several forms of the spec- 
tral density are considered in the literature, we choose an 
Ohmic bath, i.e., 



J(u})=e(ij)u, (24) 
where 0(w) is the Heaviside step function defined by 

' \ uj > 0, 



6H 



u < 0. 



(25) 
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For a proper heat conduction model a second bath of 
different temperature is coupled analogously to the op- 
posite side of the chain. The final form of the QME thus 
reads 



= -i[Hs,Ps}+V L {ps)+V R {ps 



(26) 



A full numerical investigation based on such a type of 
master equation may be found in [l2l fl3| . More generally 
the use of QME's of Redfield type like Eq. (fl~2j). respec- 
tively Eq. is common practise [U [2^, regardless of 
the well-known disadvantage, that these forms may vio- 
late the preservation of positivity of the reduced density 
operator. 

This violation is usually observed at the very early 
stage of the relaxation process only and may possibly be 
cured by a slippage of initial conditions [23, [24| . More- 
over - considering an equilibrium scenario (Pl — Pr = P) 
for a moment - the canonical Gibbs state 



-0H S 



Tr{e-^s} 



(27) 



is a stationary solution of the master equation (|26|) [251 ] - 
Whether the stationary solution of Eq. (f2"6"| is likewise 
well defined and necessarily positive for /3l ^ (3r re- 
quires further investigation. In view of these remarks 
we emphasize that we do not intend here to depart from 
Eq. (f26|) in favor of preservation of positivity itself, but 
because Eq. (|26[) can not be treated with the standard 
Monte Carlo wave function technique. Although ex- 
tended stochastic schemes exist for the solution of general 
QME's as the Redfield equation [26j], these methods turn 
out to be less efficient, in general, than the standard ap- 
proach for a Lindblad QME (see Section HTT)l . Therefore, 
it would be highly desireable to derive a Lindblad type of 
master equation for the heat conduction scenario given 
in Fig. HJ 



III. COMPLETE POSITIVITY VS. 
NONEQUILIBRIUM STEADY STATES 

In an axiomatic approach to the theory of quantum dy- 
namical semigroups the generator C of the master equa- 
tion ps = Cps can be shown to have the most general 
form [ll [3 



£ps = -i[#s,ps] 



fc 



Uk\L k p s L\ 



[L k L k , ps] 



(28) 

where [...,...]+ denotes the anti-commutator and the 
L k are operators acting on the Hilbert space of the open 
system, and the a k are non-negative numbers. A gener- 
ator of this form is said to be a Lindblad generator. It 
guarantees the preservation of the trace and of the pos- 
itivity of the density matrix. Moreover, the dynamical 
semigroup generated by C obeys the property of com- 
plete positivity [1, [27j ■ 



In general, the Redfield master equation (fT2|) . respec- 
tively (|26[) . is in Lindblad form only if further approx- 
imations are carried out. On the one hand, in many 
applications the violation of complete positivity is rea- 
sonably small and might not even show up for a large 
class of initial conditions (23[. On the other hand, it 
is always advantageous to have a Lindblad equation at 
hand, especially for numerical purposes. In fact, the 
dynamics of any QME with a generator of the form of 
Eq. (|2"5|) can be represented through a stochastic pro- 
cess in Hilbert space. This representation gives rise to 
the Monte Carlo wave function or quantum trajectory 
method 0, [ijj , which constitutes a promising technique 
to cope with open quantum systems having many degrees 
of freedom. In the present Section we consider two differ- 
ent types of derivations that aim at a QME of Lindblad 
form ((28)) for our model of heat conduction (see Fig. |T|) . 



A. Secular Approximation 

The secular approximation can be conveniently stud- 
ied by transforming to the interaction picture. In this 
representation the Redfield master equation reads 



dt 



v(pI) 



(29) 



with the Dissipator 

V(Ps)=J2 dT r ik (t, (3) [X k (t - r) pi, Xi (<)] + H.c. 

(30) 

In order to carry out the secular approximation we 
choose a suitable form for the operators X k (t) in the in- 
teraction picture. Thus - following [9( - we make use of 
a decomposition of these operators into eigenoperators 
of -Hs. With the help of projection operators n(e) that 
project onto the eigenspaces associated with the eigen- 
values e of Hs we define 



X k (w) 



e — e'—u) 



n(e)! fe n( e ') , 



(31) 



where lu represents a certain energy difference (Bohr fre- 
quency) of the system. Summing over all possible fre- 
quencies yields the desired eigenoperator decomposition 



X k =J2 x h{u) 



The operators (fBTj) obey the relation 

X\{u)) = X k (-Lu) . 



(32) 



(33) 



The interaction picture operators are now easily obtained 
from the relations 



e iAst X k (cj)e- iH ^ = e-^AfeM, 



\u)t -y-t 



Xl{u) . 



(34) 
(35) 



5 



In the case of a single system operator [see Eq. (flTl) ] 
the indices k, I in Eq. (|30[) may be dropped. Inserting 
Eqs. 421]), dSl) and (05]) into Eq. (J30j) and making use of 
the formula 

f°° P 

/ dre 1 ^-^ = nS{uj-n)+i -, (36) 

Jo lu -il 

where P indicates the Cauchy principal value. The dis- 
sipator takes the form 

^(p I s )-7r^e i (-'-)*r L (c)x 

uj.uj' 

x[X h (u,)pl,xl(uj')]+H. C . (37) 

The imaginary part of (|36|) merely leads to a small en- 
ergy shift and will silently be incorporated in the Hamil- 
tonian |T]) henceforth. By performing the secular ap- 
proximation one neglects the rapidly oscillating terms in 
Eq. (|3"T|) . Thus all terms with u>' ^ u> will be neglected 

UJ 

+ [X l (u>)pI,xI(l,)?) (38) 

Since T^,(lo) is a positive quantity for all to, Eq. (|38|) is 
of Lindblad form (|28|) as can be readily seen by a slight 
rearrangement of terms 

^ = 2nJ2^)(x L (u;) P lxl(u 1 ) 

UJ 

-±[Xt(u;)X L (w),p%} + ) . (39) 

For a heat conduction model as shown in Fig. [T] a second 
dissipator X>r(/5s) will be added to Eq. (|29|) and we end 
up with an equation analogous to Eq. (126|) . 

For systems with a non-degenerate spectrum like the 
weakly coupled spin-chains treated in this work, the SA 
decouples the equations of motion for the populations 
and coherences, i.e., for the diagonal and for the off- 
diagonal elements of /5| in the eigenbasis of the system 
Hamiltonian H$ @, H3, HE HI] , even and especially in the 
case of two baths with different temperatures. 

This decoupling gives rise to two independent systems 
of linear differential equations, one for the diagonal and 
one for the off-diagonal elements of ps in the energy rep- 
resentation. If we further assume that the stationary 
solution of Eq. (|3"5|) is unique, then this implies that the 
off-diagonal elements of ps decay and the stationary state 
/5g' is diagonal in the energy representation. 

The stationary state can therefore be decomposed as 

pf = it p n \n)(n\, (40) 

71=1 

H s \n) = e n \n) (41) 



where d denotes the dimension of the Hilbert space under 
consideration and \n) is the eigenvector corresponding to 
the eigenvalue e„ of Ha,. The expectation value of an 
operator J in the stationary state then takes the form 

(jr=Tr{p$J} 
d 

= Y,Pn (n\j\n) . (42) 

n=l 

We now require J to be an energy current operator. A 
plausible requirement on such an operator is 

(n\J\n) = 0, (43) 

since the energy eigenstates \n) must not carry energy 
currents, when the system is chain-like with open bound- 
ary conditions. It can be shown in a straight forward 
manner that our particular choice of the current opera- 
tor ((7]) satisfies the condition in Eq. (|43|) . Equation (|43|) 
enters Eq. (|42| and leads to our final result 

(J) st = . (44) 

Summarizing, the QME in the SA [cf. Eq. J39|)] fea- 
tures a steady state solution pf that is diagonal in the en- 
ergy representation, and by any plausible choice of a cur- 
rent operator, stationary currents will vanish by virtue 
of Eq. (|44|) . This prediction is most unphysical in view 
of the nonequilibrium scenario that is depicted in Fig.[T]. 

In this sense the QME in the SA proves inappropriate 
for our purposes. This failure of the SA is physically plau- 
sible since, as already mentioned, it only applies if the 
differences of the Bohr frequencies are large compared to 
the relaxation rates of the system. In chainlike systems, 
however, differences of Bohr frequencies between energy 
levels within the bands will go to zero with chainlengths 
going to infinity (even in the limit of strong couplings). 
Thus, the differences of the Bohr frequencies will eventu- 
ally be too small to justify the application of the SA in 
a nonequilibrium scenario. 

B. Weak Internal Coupling Approach 

In this Section we shall derive a Lindblad QME that 
retains the property of Eq. ()26|) of featuring a nonequi- 
librium steady state with non- vanishing energy current, 
and is thus suitable for the description of energy trans- 
fer. To this end, we assume that for weak couplings of 
neighboring spins X(—t) [cf. Eq. JT5J)] is approximately 
given by the free dynamics (see [301), i- e -> 

X{-t) = ( :// WV" ' ~ e-'^'^Xe 16 ^ . (45) 

We would like to stress that the weak coupling assump- 
tion must be considered fundamental for the very model 
by virtue of Eqs. ([7]) and (|2l)|) . The time dependent ver- 
sion of the operator (| 1 T[) thus takes the simple form 

X h (-r) ~ e- inT a ( + } + e inT . (46) 
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This approximation enters Eq. (|T8|) and yields 



drr(r,/3 L )x 



x e 



,-if2Tr-(l) ~ 



iQr r-(l) 



r(l)' 



(47) 
H.c. 



with r(r, given by Eq. (j2Tj) . By making use of 
Eq. (|36p and neglecting the imaginary part (see remarks 
in Sec. IIII A[) we are left with 



(48) 



To proceed we introduce a complete set of operators 
Fi acting on the Hilbert space of the first spin. It is con- 
venient to choose this operator basis to be orthonormal 
with respect to the Hilbert-Schmidt scalar product 



(F i ,F j )=Tr{F}F j } = S i: 



(49) 



A suitable set of operators satisfying these conditions is 
given by the set {cr+, <5\_, <r z /2, 1/2}. With the help of 
these operators Eq. (|48|) can be cast into the form 



2 ! 



k,l=l 



(50) 



with 



ft 

Pi — <J + 



Fi = a 



1(2) 
1(2) 



l(n) 



(51) 



The coefficient matrix 7 = (7/0;) introduced in Eq. (|50|) 
reads 

/ 2r L (^) r L (Q) + r L (-n)\ 
7 -Mr L (fi) + r L (-n) 2r L (-fi) )• [bZ) 



A given dissipator of the form of Eq. (|50[) is in Lindblad 
form if and only if the coefficient matrix 7 is positive (see, 
e.g., [§]). However, the matrix given by Eq. |52|) is not 
positive because 711 > and 



det( 7 ) = -7r[r L (fi) - r L (-r!)] 2 < o 



(53) 



since r^ifi) £ R. Hence, we conclude that the dissipa- 
tor derived above is not in Lindblad form. 

We suggest the following strategy to bring the dissipa- 
tor given by Eq. (|50[) into Lindblad form by only mini- 
mal modifications, without invoking the SA. The idea is 
to separate off the largest possible contribution to 7 that 
leads to a generator in Lindblad form. To this end, we 
first decompose the coefficient matrix as 



where 



7 = 1A + IB 



In M 
~"' = ' 7 I M 722 

712 - M 
721 - M 



7B = 7T 



(54) 

(55) 
(56) 



Of course, such a decomposition holds true for any value 
of the parameter M. We fix this parameter by the re- 
quirement that the determinant of the matrix 7^4 van- 
ishes, 



det ( 7A ) = 0. 



(57) 



This condition implies that is positive, having one 
zero eigenvalue. Obviously, the condition (|57|) is satisfied 
if we choose 



m = 2Vr L (tt)r L (-ft) . (58) 

With this choice the matrix 75 becomes 



where 



7s = ttk J(tt) I ° q 



A = N+ - - VA^T~/V 



(59) 



(60) 



with the Planck distribution N = N(J3,Q.) [cf. Eq. (|23"])]. 
We see that the non-zero elements of 75 rapidly vanish 
as the temperature and, likewise, N goes to infinity 



lim ,4 = 0. 



(61) 



For not too low temperatures the quantity A is small, and 
we may neglect the contribution from jb to the dissipa- 
tor, ending up, finally, with a master equation in Lind- 
blad form. Neglecting jb in an appropriate parameter 
range can be considered as a minor invasion in compari- 
son to the numerous assumptions that lead to Eq. (|50[) . 
Summarizing, the final form of the QME in the weak 
internal coupling limit thus reads 

^ = -i{H Sl p s }+V L (p s ) + V R (p s ), (62) 
where the dissipator for the left heat bath is given by 



Pl(/3s)= £(7a)«( 
fe,i=i 



F k p s F}--[F}F k ,p s ] + ), (63) 



and the dissipator 2?r,(/3s) for the right heat bath is de- 
fined correspondingly. 

Figure [2] demonstrates the excellent match between 
the predictions of the Redfield master equation (TT!?| and 
those of the Lindblad form (f62|) that has been derived in 
this Section. Furthermore the results from a stochastic 
wave function simulation of Eq. (|62|) are depicted. We 
refer the reader to 0, [ItJ for details of the method. In- 
vestigations on larger chain models now become feasible 
and more detailed numerical results will be presented in 
a forthcoming paper. 

Thus, the Eq. (|6"2"]) contains both vital properties: On 
the one hand it is in Lindblad form, preserving all prop- 
erties of the density matrix and being, furthermore, well 
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FIG. 2: Energy current (J*- 1,2 - 1 ) in a chain of n = 3 weakly 
coupled spins, driven by heat baths of different temperature 
{/3 t = 0.41, /3 R = 1.39, X — k — 0.01, n = 1). Solid line: 
computed on the basis of the Redfield master equation (|26|l by 
exact diagonalization in Liouville Space. Filled circles: weak 
internal coupling approximation (|62[) (again exact diagonal- 
ization) . Circles: Numerical solution of Eq. (|62|l by means of 
a stochastic wave function simulation (10 5 realizations). 

suited for all types of stochastical methods. On the other 
hand it features a nonequilibrium steady state with an 
energy current flowing through the system. 

The neglect of the interaction between spins in Eq. ([43]) 
leads to the effect that the damping concerns the hypo- 
thetically uncoupled system only. All effects that arise 
from different internal coupling schemes will only show 
up in the coherent part of the QME. However, the dy- 
namics will not differ in the type of damping. One possi- 
ble way to reintroduce the effects of the internal coupling 
might be a perturbational approach [3(| , but again posi- 
tivity violation will arise most probably without further 
modifications. Moreover we would like to point out that 
the neglect of the off-diagonal entries in Eq. ([52]) directly 
yields a Lindblad QME which is of the same type that 
has previously served for investigations on heat trans- 
port in small quantum systems [14j but hitherto lacked 
a microscopic derivation. 

IV. CONCLUSIONS 

The purpose of this article was, on the one hand, to 
present a model of an open system for which the stan- 



dard microscopic derivation of a Lindblad quantum mas- 
ter equation (QME) fails, in the sense that it is inappro- 
priate to model nonequilibrium steady states when the 
system is subject to a temperature gradient induced by 
baths of different temperature. 

On the other hand, we have proposed an alternative 
derivation that does yield a Lindblad QME and can prop- 
erly describe energy currents in the stationary state of 
weakly coupled spin-i Heisenberg chains. The usefulness 
of a QME of Lindblad-form can not be overemphasized, 
because it is the key link to efficient numerical investi- 
gations by means of standard stochastic wave function 
methods. 

Stress was laid on the fact that the secular approxima- 
tion (SA) performed on a quantum master equation for 
a non-degenerate, one-dimensional system in a split bath 
scenario with temperature gradient (see Fig. []} leads to a 
stationary state which is diagonal in the energy eigenba- 
sis of the system. Therefore, no energy current is flowing 
from the hot towards the cold bath, a situation which 
is highly unphysical. Hence, the use of a QME in the 
SA is inappropriate for the description of this type of 
nonequilibrium scenarios. 

An open question remains whether the secular approx- 
imation generally wipes out all local coupling aspects as 
considered in this paper. Our studies have hitherto been 
restricted to one-dimensional spin chains only. One pri- 
mary future goal therefore is to carry out further system- 
atic investigations on the question whether our findings 
are extendable to larger classes of open systems and cou- 
pling models. 
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